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There is growing experimental interest in coupling cavity pho¬ 
tons to the cyclotron resonance excitations of electron liquids in 
high-mobility semiconductor quantum wells or graphene sheets. 
These media offer unique platforms to carry out fundamental stud¬ 
ies of exciton-polariton condensation and cavity quantum electro¬ 
dynamics in a regime in which electron-electron interactions are 
expected to play a pivotal role. Focusing on graphene, we present a 
theoretical study of the impact of electron-electron interactions on 
a quantum Hall polariton fluid, that is a fluid of magneto-excitons 
resonantly coupled to cavity photons. We show that electron- 
electron interactions are responsible for an instability of graphene 
integer quantum Hall polariton fluids towards a modulated phase. 
We demonstrate that this phase can be detected by measuring the 
collective excitation spectra, which soften at a characteristic wave 
vector of the order of the inverse magnetic length. 

Fluids of exciton polaritons, composite particles resulting from coupling 
between electron-hole pairs (excitons) in semiconductors and cavity photons, 
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have been intensively investigated over the past decade Because 

of the light mass of these quasiparticles, exciton polariton fluids display 
macroscopic quantum effects at standard cryogenic temperatures BEIE], 
in stark contrast to ultracold atomic gases. Starting from the discovery of 
Bose-Einstein condensation of exciton polaritons in 2006 [1], these fluids have 
been the subject of a large number of interesting experimental studies explor¬ 
ing, among other phenomena, superfluidity PE], hydrodynamic effects |7], 
Dirac cones in honeycomb lattices [8], and logic circuits with minimal dissi¬ 
pation [9]. 

The isolation of graphene nni-a two-dimensional (2D) honeycomb crys¬ 
tal of carbon atoms— and other 2D atomic crystals m including transition 
metal dichalcogenides [mE] and black phosphorus HI, provides us with an 
enormously rich and tunable platform to study light-matter interactions and 
excitonic effects in 2D semimetals and semiconductors. Light-matter interac¬ 
tions in graphene in particular have been extensively explored over the past 
decade with both fundamental and applied motivations dSl DEI HZl EHl Ej. 
Recent experimental advances have made it possible to monolithically inte¬ 
grate graphene with optical microcavities [T9| EQ], paving the way for fun¬ 
damental studies of cavity quantum electrodynamics at the nanometer scale 
with graphene as the active medium. Progress has also been made using 
an alternate approach applied previously to conventional parabolic-band 2D 
electron liquids in semiconductor quantum wells |2T] by coupling graphene 
excitations to the photonic modes of a Terahertz (THz) metamaterial formed 
by an array of split-ring resonators |22) . 

When an external magnetic field is applied to a 2D electron liquid in a 
GaAs quantum well j23] or a graphene sheet 1^ . and electron-electron 
interactions are ignored, transitions between states in full and empty Landau 
levels (LLs) are dispersionless, mimicking the case of atomic transitions in a 
gas. The cyclotron resonance of a 2D quantum Hall fluid can be tuned to 
resonance with the photonic modes of a cavity or a THz metamaterial j21) . 
thereby establishing the requirements for cavity quantum Hall electrodynam¬ 
ics (cQHED). Cavity photons have already been used to carry out spectro¬ 
scopic investigations of fractional quantum Hall fluids |26]. cQHED phe¬ 
nomena present several important twists on ideas from ordinary atom-based 
cavity quantum electrodynamics because in this case interactions between 
medium excitations are strong and long-ranged. Furthermore the active 
medium can be engineered in interesting ways, for example by using, in¬ 
stead of a single 2D crystal, van der Waals heterostructures I271I2HIE9] or 
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vertical heterostructures which include both graphene sheets and ordinary 
semiconductor quantum wells [301 ED. 

In this Article we show that electron-electron (e-e) interactions play a ma¬ 
jor qualitative role in graphene-based cQHED. Before describing the technical 
details of our calculations, let us briefly summarize the logic of our approach. 
The complex many-particle system of electrons in a magnetic field, interact¬ 
ing between themselves and with cavity photons, is treated within two main 
approximations. We use a quasi-equilibrium approach based on a micro¬ 
scopic grand-canonical Hamiltonian and treat interactions at the mean-field 
level. We critically comment on these approximations after the Results sec¬ 
tion. Our approach is similar to that used in Refs. [321 E3l IMl, except that 
simplifications associated with Landau level quantization allow more steps 
in the calculation to be performed analytically. 

The problem of finding the most energetically favorable state of a graphene 
integer quantum Hall polariton fluid (QHPF) is approached in a variational 
manner, by exploiting a factorized many-particle wave function. The latter 
is written as a direct product of a photon coherent state and a Bardeen- 
Cooper-Schrieffer state of electron-hole pairs belonging to two adjacent LLs. 
We find that e-e interactions are responsible for an instability of the uniform 
exciton-polariton condensate state towards a weakly-modulated condensed 
state, which can be probed experimentally by using light scattering. We 
therefore calculate the collective excitation spectrum of the graphene QHPF 
by employing the time-dependent Hartree-Fock approximation. We demon¬ 
strate that the tendency to modulation driven by e-e interactions reflects into 
the softening of a collective mode branch at a characteristic wave vector of 
the order of the inverse magnetic length. 

Results 

Effective model. We consider a graphene sheet in the presence of a 
strong perpendicular magnetic field B = Rz [351136] . We work in the Landau 
gauge with vector potential Aq = —Byk. The magnetic field quantizes 
the massless Dirac fermion (MDF) linear dispersion into a stack of LLs, 
^\,n = Xhwc\/n, which are labeled by a band index A = ± which distinguishes 
conduction and valence band states and an integer n = 0,1, 2,.... Here 
ojc = is the MDF cyclotron frequency [35l[36], — c/300 the Dirac 

band velocity (c being the speed of light in vacuum), and = \/hc/(eR) ~ 
25 nm/^R [Tesla] is the magnetic length. The spectrum is particle-hole 
symmetric, i.e. = —c+,n for each n. Each LL has macroscopic degeneracy 
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N = N^S /where iVf = 4 is the spin-valley degeneracy and 
S' = is the sample area. 

In this Article we address the case of integer hlling factors, which we 
expect to be most accessible experimentally. Because of particle-hole sym¬ 
metry, we can assume without loss of generality that the chemical potential 
lies in the conduction band between the n = M and n = M -I- 1 Landau 
levels. When the energy hw of cavity photons is nearly equal to the cy¬ 
clotron transition energy VLm = £+,m+i — ^+,m = hudVM + 1 — \/M), the 
full fermionic Hilbert space can be effectively reduced to the conduction-band 
doublet M, M + 1. 

We introduce the following effective grand-canonical Hamiltonian: 


— '^ph + '^mat + 'Hint “ he^e — hX^X • (1) 


The hrst term, "Hph, is the photon Hamiltonian, "Hph = J2qu 
where j, {ciq,u) creates (annihilates) a cavity photon with wave vector q, 
circular polarization n = L, R, and frequency ujq = being 

the cavity dielectric constant and c the speed of light in vacuum. 

The second term in Eq. ([^, 'Hma.t, is the matter Hamiltonian, which 
describes the 2D MDF quantum Hall fluid, and contains a term due to e- 
e interactions. This Hamiltonian is carefully derived in the Supplementary 
Note 1. In brief, one starts from the full microscopic Hamiltonian of a 2D 
MDF quantum Hall fluid |36j . written in terms of electronic held operators 
c\^n,k,i- Here, A = ± is a conduction/valence band index, n is a LL index, 
k = —L/{21\) + {27r/L)j with j = is the eigenvalue of the x- 

direction magnetic translation operator, and is a fourfold index, which 
refers to valley {K, K') and spin O', J,) indices. All the terms that involve 
held operators c\^n,k,i, acting only on the conduction-band doublet 

M, M + 1 are then treated in an exact fashion, while all other terms are 
treated at leading order in the e-e interaction strength ra- 

The third term. Hint, describes interactions between electrons and cavity 
photons, which we treat in the rotating wave approximation. This means that 
in deriving Hint we retain only terms that conserve the sum of the number 
of photons and the number of matter excitations. Details can be found in 
the Supplementary Note 2. It is parametrized by the following light-matter 
coupling parameter 




hjjJr 


2KrL.hUn 


( 2 ) 
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In Eq. ^ Lz L is the length of the cavity in the z direction (V = 
LzLF' is the volume of the cavity). In what follows we consider a half¬ 
wavelength cavity, setting u = nc/{Lzy/K~r)- Consequently [25], g = go = 
hujc\/aQED/, where agso = e^/{hc) ~ 1/137 is the QED hne- 
structure constant. 

Finally, in Eq. Q we have introduced two Lagrange multipliers, /ig and 
/ix, to enforce conservation of the average number of electrons and excita¬ 
tions [2B]- Ne is the electron number operator in the M, M -f 1 reduced 
Hilbert space, while Nx = iVph -|- is the operator for the number of mat¬ 
ter excitations (excitons). The value of the chemical potential fig should be 
hxed to enforce {'ijj\Ne\'ip) = M. At zero temperature this condition is simply 
enforced in the variational wave function dehned below. 

Variational wave function and spin-chain mapping. To hnd the 
ground state of the Hamiltonian Q, we employ a variational approach in 
which the many-particle wave function \'ip) is written as [3^ 133] a direct 
product of a photon coherent state and a Bardeen-Cooper-Schrieffer state of 
electron-hole pairs belonging to the M, M + 1 conduction-band doublet: 


IV’) = 


X 


exp exp {yMa oJj l) 

JJ Jcos(6'fc/2) -F sin{ek/2)c\ j^z,i^^^^c+^M,k,i 


IVo) , 


( 3 ) 


where |'0o) is the state with no photons and with the M-th LL fully occu¬ 
pied. In writing Eq. ([^, we have allowed for phase-coherent superposition 
of electron-hole pairs with fc-dependent phases 0^ and excitation amplitudes 
sin(6'fc/2), to allow for the emergence of modulated QHPF phases driven by 
e-e interactions. Eq. ([^ can be written in terms of polariton operators, as 
shown in Supplementary Note 3. The variational parameters {(^k}, 

and a can be found by minimizing the ground-state energy . We 

introduce the following regularized energy (per electron) |40j 


S = S{{(j)k},{0k},a) 


i'lpl'Hl'ip) - (V’olTfiV’o) 
Af 


( 4 ) 


The variational wave function ([^ and the functional T({0fc}, {0^}, a) can 
be conveniently expressed in terms of the /c-dependent Bloch pseudospin 
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orientations: 


n(/c) = = [sin(6'fc)cos((?!)fc),sin(6'fc)sin(0fc),-cos(6'fc)]'^ , 

iVf ^ 

(5) 

where ak,^ = Y.n,n'=M,M+i 'Tnn'C^+^n,k,^^+y,k,i and r = (ti, Ta, rg) is a 3D vector 
of Panli matrices acting on the M, M + 1 doublet. The variational wave 
function then becomes 

IV’) = exp exp oIj l) 

X exp z '^9km{k) ■ (Tk,^/2^ l-^o) , (6) 

where m(/c) = (sin(0fc), — cos(0fc),0) is a unit vector orthogonal to n(/c) 
and I'^o) contains all pseudospins oriented along the —z direction. Since 
exp [—z6*A:ni(fc) ■ <Jk,^/‘2] acts as a rotation by an angle 9^ around m(/c), we can 
interpret the matter part of jV) as a state in which every pseudospin labeled 
by (k,^) is rotated accordingly. The unit vector n{k) in Eq. (|^ denotes the 
hnal pseudospin direction at each k = 1. . The string {n(fc)}fc of 

unit vectors can be viewed as a set of classical spins on a one-dimensional 
chain whose sites are labeled by the discrete index fc, as in Fig. 

In the same notation, 


T = e — eo -l- hwnx 


(7) 


where 

1 1 1 f ^ 

e = -g-^^B-n{k) +--^^l^Je{k-k')neik)ne{k') 

-I- T>{k — k') ■ [n(/c) X n(fc')] j" , ( 8 ) 

fuD = ha;q=o — /zx = ^ — hx, and 


—AT- 



+ 


1 n2,{k) -|- 1 


(9) 
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In Eq. ([^, the quantity hw plays the role of a Lagrange multiplier, eo is 
a reference energy which is defined so that that £ = Q when nx = 0, and 
££ = (j7i, d/ 2 , jTs) and = (0,1^2,0) are the symmetric and antisymmet¬ 
ric, i.e. Dzyaloshinsky-Moriya (DM) H2] . interactions between Bloch 
pseudospins. Explicit expression for X>, dT”, and eo are provided in Sup¬ 
plementary Note 4, together with plots of the Fourier transforms J^{q) = 

and:D 2 (g) = Zlfc ^ 2 (fc) exp in Sup¬ 

plementary Figure 2. In Eq. 0 we note a photon contribution riph = 
('0|-Aph|'0)/A/' = |ap and an exciton contribution, nex = ('0|A"ex|'0)/A^ = 
A/"^^ sin^(6'fc/2). It is somewhat surprising that DM interactions appear 

in our energy functional ([^, since these require spin-orbit interactions and 
appear when inversion symmetry is broken. Our microscopic Hamiltonian 
([^ does not contain neither SOIs nor breaks inversion symmetry. In the 
next Section we discuss the origin of pseudospin DM interactions. 

Each of the terms in the expression ([^ for e = e({0fc}, {0^}, a) has a clear 
physical interpretation. The first term on the right-hand side is the energy 
of a set of independent ID Bloch pseudospins in an effective magnetic field 
with the usual Rabi coupling and detuning contributions: 


B 


-\/2Re(a). V2Im(a), (A - aee)/( 25 ) , 


( 10 ) 


where A = hw — {Qm + Age) is the detuning energy with Agg a correc¬ 
tion due to e-e interactions between electrons in the M, M + 1 doublet and 
electrons in remote occupied LLs na El sg (see Supplementary Note 4). 
Because the MDF model applies over a large but finite energy interval, we 
need to introduce an ultraviolet cut-off Umax on the LL labels n of occupied 
states with A = — 1. Our choice for Umax is explained in the Supplementary 
Note 5. It is easy to demonstrate that Agg depends logarithmically on U ma x: 
Agg = (aggO^/S) [In(nmax) + Cm\ where Ogg = e^/{Krhvji) is the graphene 
fine structure constant jl6] and Cm is an ultraviolet-finite constant. For 
M = 1 we find that Ci ~ —2.510, in agreement with earlier work [IS]- The 
correction Agg to the cyclotron transition energy is related to the extensively 
studied naiiHiiin] renormalization of the Dirac velocity ud due to exchange 
interactions which also occurs in the absence of a magnetic field. The quan¬ 
tity Ogg involves only e-e interactions within the M, M + 1 doublet (see Sup¬ 
plementary Note 4). For M = 1 we find Ogg = agghcUg(63/32 -|- y/2)y/Ti/‘i2. 
The second term in Eq. (|^ describes interactions between Bloch pseudospins, 
which originate microscopically from matter-coherence dependence in the e-e 
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interaction energy. At long wavelength these interactions stiffen the polari- 
ton condensate collective mode dispersion and support superffuidity. In the 
absence of a magnetic field their role at shorter wavelengths is masked by 
increasing exciton kinetic energy |50| . 

Pseudospin DM interactions. In the QHPF exciton fluid kinetic en¬ 
ergy is quenched and, as we explain below, DM exciton-exciton interactions 
play an essential role in the physics. We therefore need to understand why X> 
is finite. We start by observing (see Supplementary Note 4) that X> contains 
direct D 2 ,d and exchange 1 ^ 2 ,x contributions, which a) are of the same order 
of magnitude and b) have the same sign. We can therefore focus on the direct 
contribution, which has a simple physical interpretation as the electrostatic 
interaction between two charge distributions that are uniform along the x 
direction and vary along the y direction, i.e. 

V 2 ,d{k - k')n^{k)n^{k') = J dydy'ln Pz{y,k)p^{y\k') 

( 11 ) 

where 

Pz{y.k) = -0^(2/-A:4) 

and 

Px{y,k) = -enx{k)[w+^M4>M+i{y - ki%)(j)M{y - kis) 

+ w_^M(j)M{y - ki%)(j)M-i{y - ki%)] . (13) 

Here 4>n{y) with n = 0,1,2,... are normalized eigenfunctions of a one¬ 
dimensional harmonic oscillator with frequency ojc and = y/l ± 5np 

captures the property that the pseudospinor corresponding to the n = 0 LL 
has weight only on one sublattice |36j . 

We now use a multipole expansion argument to explain why T> 2 ^d{k — k') ^ 
0. We first note that pz{y,k) has zero electrical monopole and dipole mo¬ 
ments but finite quadrupole moment Q{k) = —ei%QMnz{k) = — 

dwifl/‘^)nz{k). On the other hand, px{y,k) has zero electrical monopole but 
finite dipole moment d{k) = —eiBdMnx{k) = + l)/2 -|- 

\/^^dj2]nx{k) . Using a multipole expansion, it follows that the leading con¬ 
tribution to Eq. ([TT| is the electrostatic interaction between a line of dipole 
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moments extended along the x direction and centered at one guiding center 
and and a line of quadrupole moments centered on the other guiding center. 
It follows that 


V2Ak - k') 


2e^ QndM 
eAfLil {k - k'f ■ 


(14) 


The interactions are antisymmetric, he. their sign depends on whether the 
dipole is to the right or to the left of the quadrupole. The direct contribu¬ 
tions between like pseudospin components which contribute to JT” are sym¬ 
metric because they are interactions between quadrupoles and quadrupoles 
or dipoles and dipoles. 

Alert readers will have noted that only the y-direction DM interaction is 
non-zero, T>{k — k') = 'D 2 {k — k')y . In contrast, the usual DM interaction |41l 
142] is invariant under simultaneous rotation of orbital and spin degrees of 
freedom. This is not the case for pseudospin DM interactions: the property 
that only the y component of is non-zero can be traced to the property 
that, for a given sign of pseudospin n^{k), the charge distribution Px{y,k) 
in Eq. (13) changes sign under inversion around the guiding center (i.e. y —)■ 
-y + 2kil). 

Linear stability analysis of the uniform fluid state. We hrst assume 
that the energy functional is minimized when 9k and 0^ in Eq. ^ are k- 
independent, i.e. 9k = 9 and (j)k = 4> for every k. The functional £ then 
simplifles to 


£{(/), 9, a) = -f (ho; — A) sin^(6'/2)-|-Oee sin"‘(6*/2) 

+ 2^2^71 a\ cos(0 -I- arg(a)) sin(6'/2) cos(6'/2) . 


(15) 


The hrst term on the right hand side of Eq. (15), which is proportional to 


|af, is the free photon energy measured from the chemical potential px- The 
second term, which is proportional to sin^(6'/2), is the free exciton energy 
(as renormalized by e-e interactions, which enter in the dehnition of A). 
The third term, which is proportional to sin^(6'/2), is the exciton-exciton 
interaction term. Finally, the term in the second line, which is proportional 
to the Rabi coupling \/2g, describes exciton-photon interactions. 

We seek for a solution of the variational problem 6£ = 0 characterized 
by non-zero exciton and photon densities. For this to happen, the common 
chemical potential px needs to satisfy the following inequality: 


\huj{huj — A -I- Oee)! < 2g^ + hwae 


(16) 
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When this condition is satisfied, the solution oi 5£ = Q is given by 


arg(a) = n - <t> , 


( 17 ) 


lap = 




1 - 


hu{hu — A + Oee) 
2g‘^ + huttee 


and 


cos(6*) = 


hu{hu — A + Oee) 


(18) 


(19) 


{2g^ + huttee) 

The common chemical potential px must be adjusted to satisfy ux = [1 — 
cos(6')]/2 + |ap. In the spin-chain language introduced above this state is 
a collinear ferromagnet in which all the classical spins {n(/c)}fc are oriented 
along the same direction, as in Fig. [^a). Note that, as expected, the energy 
minimization problem does not determine the overall phase of the condensate. 

We now carry out a local stability analysis to understand what is the 
region of parameter space in which this polariton state is a local energy 
minimum. A minimum of S, subject to the constraint on the average density 
nx of excitations, is also a minimum of the functional e({0fc}, {Ok}, «) defined 
in Eq. ([^ with |a| not considered as an independent variable but rather 
viewed as a function of the variational parameters {6^} through the use of 

Eq. (9), i.e. with |a| -)■ |a„x({6'fc})| = ^rix - With 

this replacement, 

e({0fc},{0fc},a)||„|^l„(g^)l (20) 

becomes a functional of 2Af^ + 1 independent variational parameters, which 
can be arranged, for the sake of simplicity, into a vector w with components 
w = (arg(a),0i,...,0Ar^,0i,...,0ArJ^. 

In this notation, the extremum discussed above can be represented by 
the vector wq = (tt — 0,6*,..., 6*, 0,..., 0)"^. We have checked that wq is a 
solution of the equation Vw?^(w) = 0. Whether wq is a local minimum or 
maximum depends on the spectrum of the Hessian 


-Amn(wo) = 


02 


e w 


dWmdWr, 


( 21 ) 


w=wo 


which is a {2N'^ 4-1) x (2A/(i -|- 1) symmetric matrix. 

The homogeneous polariton fluid phase is stable only if iFmn(wo) has 
no negative eigenvalues. The stability analysis is simplified by exploiting 
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translational symmetry to classify state flnctnations by momentnm. Stability 
phase diagrams for M = 1 and M = 2, constrncted by applying this criterion, 
are plotted in Fig. [^for two different valnes of the cavity dielectric constant 
Kj-. In this fignre, white (grey-shaded) regions represent the valnes of the 
detnning A and density nx of total excitations for which the homogeneons 
flnid phase is stable (nnstable). As expected, by increasing (i.e. redncing 
the importance of e-e interactions) the stable regions expand at the expense 
of the nnstable ones. Note that the instability displays an intrigning reentrant 
character and that it can occnr also when matter and light have comparable 
weigth, i.e. when riex ~ "n-ph- 

We have checked that the root of instability of the homogeneons flnid 
phase is e-e interactions. More precisely, it is possible to see that in the 
absence of DM interactions— i.e. when = 0 in Eq. ([^— the instability 
disappears. Symmetric interactions, however, still play an important qnan- 
titative role in the phase diagrams, as explained in Snpplementary Note 6. 
The physics of these phase diagrams is discnssed farther below where we 
identify the phase diagram bonndary with the appearance of soft-modes in 
the nniform polariton flnid collective mode spectrnm. Stable phases occnr 
only if a) > 0 (i.e. px < hw). We remind the reader that this condition on 


Hx is additional to the one given in Eq. (16) above. 


Elementary excitations of the polariton fluid. We evalnate the el¬ 
ementary excitations of the nniform polariton flnid by linearizing 

the Heisenberg conpled eqnations of motion of the matter Bloch psendospin 
and photon operators nsing a Hartree-Fock factorization for the e-e interac¬ 
tion term in the Hamiltonian (see Snpplementary Note 7). The collective 
excitation energies diagonalize the matrix 


/ 

hcUq 

0 



\ 


0 

-hAq 

9x' 






(q) 

— z[DmF + ‘^Jy'y 

(q)] 

V 

—2igx\q 

-2w:-,q 

z[DmF + 2Ja;/a;/(q)] 

i2Jyix' (q) 

/ 


M = 


( 22 ) 

The hrst two-components of eigenvectors of M correspond to photon cre¬ 
ation and annihilation, and the third and fonrth to rotations of the Bloch 
psendospin in a plane (denoted by x'-y' in the Snpplementary Note 7) or¬ 
thogonal to its gronnd state orientation. 


In Eq. (22) 


hcjq = 




11 






(23) 


gx',q = e^'^cos{9)gg/V2, gy>^g = ie*%/A/2, 


n 


MF — 


hw 


2J-i(0) , 


Jx'x'{oi) = [j7i(g)sin^(v3q-0) + J^2(?)cos^(</3q-0)]cos^(6') 


+ J3(?)sin2(0) , ^24) 

4'y'(q) = J 2 {.q) sin^(v5q - 0) + Ji{,q) cos^(</3q - 0) , (25) 

sRe[Ja;v(q)] = sin(2(^q - 20) cos(6')[^2(g) - ^i(g)]/2 , (26) 

'^m[J^,yi{q)\ = V{q) cos{g)g - cf)) sin(6') , (27) 

and q = [qcos{ipg),qsm{ipg)]'^. 


The solution of the eigenvalue problem yields two hybrid modes that 
can be viewed as lower polaritons (LP) and upper polaritons (UP) that are 
dressed by the condensate and have strong mixing between photon and mat¬ 
ter degrees of freedom at isq < vd-^/k^/c -C 1. Fig. [^illustrates the disper¬ 
sion relations of these two modes for M = 1. For wavelengths comparable to 
the magnetic length, qis ~ 1, the UP mode has nearly pure photonic char¬ 
acter, while the LP mode is a nearly pure matter excitation with a dispersion 
relation that is familiar from the theory of magnon energies in systems with 
asymmetric DM exchange interactions |51j : 


Dq ^ 2^?7r (q)] 

+ y^[DmF + 2Jx'a;'(q)][f2MF + x'y'{ q^]Y ■ (28) 


Fig. ||(b) shows the LP dispersion relation for three different polar an¬ 
gles (fg. In all cases, a local roton-like minimum occurs at a wave vector 
2|^B < q < 3/£s. The global minimum of the LP dispersion occurs along 
the direction ipg = cj), where the impact of DM interactions is strongest, 
i.e. ‘^m[Jx'y'{q)] is maximum—see Eq. (27). The mode energy vanishes, and 
a Hessian eigenvalue crosses from positive to negative signalling instability, 
when 


[Dmf + 2J3;/a;/(q)] [Dmf + 2Jy'j/'(q)] < 


Jx'y 


(q)P 


(29) 


Since the LP mode becomes unstable at a finite wave vector rNj l/is, we 
conclude that the true ground state spontaneously breaks translational sym¬ 
metry. We emphasize that softening of collective modes in quantum Hall 
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fluids can be experimentally studied, for example, by inelastic light scatter¬ 
ing |52]. 

Modulated phase of QHPFs. Motivated by the properties of mag¬ 
netic systems with strong asymmetric spin interactions |5Tj, we seek broken 
translational states in which the Bloch pseudospins execute a small ampli¬ 
tude spiral around a mean orientation, as in Fig. [^b). This is a state in 
which 6k, 4>k have a rather simple ^-dependence of the form; 


6k = 6 + u cos{£%Q*k) , 

4>k = ~ vsm{i‘^^Q*k + Lf) . 


(30) 


Eq. (30) physically describes a small-amplitude spatially-periodic contribu¬ 
tion to the uniform condensate state ([^ with 6k = 6 and (pk = 0- One should 
therefore not confuse the condensed state described by Eqs. (|^ and ( |30| with 
a uniform condensate in which electrons and holes form pairs with a hnite 
center-of-mass momentum |53) . 

Because the form factors of electrons in M and M -I- 1 LLs differ, this 
state has non-uniform electron charge density with periodicity (27r/Q*)y. 
The Fourier transform of the density variation 


(5nq = 


£y [(V^|T'^(r)T(r)|^/;) - (V’o|d'^'(r)T(r)|V'o)] , (31) 


is non-zero only for q = (p),nQ*), where u is a relative integer. In Eq. (31) 
'J'W = EA,n,fc,€(r|A,n, k)cx,n,k,i with 


{r\X,n,k) = 


e*"" f W_,r^cPn-l{y - £lk) 


V >^W+^n4>n{y - £lk) 


(32) 


is a held operator that creates an electron at position r [36]. 

For this form of variational wave function we have hxed 6, p, a, u, v, Q*, 
and (p by minimizing £. A summary of our main numerical results for u, v, 6, 
and Q* is reported in Fig. |^for two values of the detuning A. Minimization 
yields (p = 0, 0 = —1^/2, and arg(Q;) = tt — 0. The dependence of |q;| on 
nx is given by |a| = — [1 — cos(6*) Jo(m)] /2, where Jo{x) is the Bessel 

function of order zero. Figs. |^a) and (b) illustrate the weak dependence of 
the characteristic wave number Q* on the density nx of total excitations. In 
Fig. [^we report, for each value of nx, the ratio 


Te 


Cm ~ Ch 
|eh - eo| 


(33) 


13 








The numerator in Eq. (33) is the difference between the energy of the con¬ 
densed modulated phase, em, described by Eq. (30), and that of the condensed 
homogeneous phase, eh, described by Eqs. (17)-(19). In the condensed mod¬ 
ulated phase, it follows that < 0. The denominator in Eq. (33) is the 
condensate energy in the homogeneous phase. In Fig. [^we clearly see that, 
depending on the detuning A and the density of excitations nx, the modu¬ 
lated phase comes with a condensate energy gain r^. in the window 5%-15%, 
with values of the photon fraction that are well above 10%. 

Discussion. In this Article we have made two major simplifying ap¬ 
proximations that deserve a detailed discussion. We have i) used a quasi¬ 
equilibrium approach based on a grand-canonical Hamiltonian and ii) treated 
electron-electron interactions at the mean-field level. 

i) Exciton-polariton condensates differ from ultracold atomic gases in that 
the condensing quasiparticles have relatively short lifetimes, mainly because 
of photon losses in the cavity or metamaterial. External optical pumping 
is therefore needed to maintain a non-equilibrium steady state. It has been 
shown lai that the resulting non-equilibrium steady state can be approxi¬ 
mated by a thermal equilibrium state when the thermalization time is shorter 
than the exciton polariton lifetime. Equilibrium approximations have been 
successfully used in the literature to describe exciton-polariton fluids in semi¬ 
conductor microcavities [21 ESI ESI ESI ES]- Experimental studies in GaAs 
quantum wells have shown that the thermalization time criterion is satis¬ 
fied above a critical pump level 153 and that polariton-polariton interactions 
(which are responsible for thermalization) are strong |58j . We assume below 
that a similar thermal equilibrium state can be achieved in graphene QHPFs. 
Because polaritons interact more strongly when they have a larger excitation 
fraction, quasi-equilibrium polariton condensates are expected to be more 
accessible experimentally when the cavity photon energy is higher than the 
bare exciton energy, i.e. at positive detuning. 

ii) The possibility of non-mean-field behavior in the matter degrees of 
freedom is an issue. Mean-field theory is accurate for dilute excitons at low 
temperatures [52], but could fail at high exciton densities. In particular, 
the modulated phase we have found may undergo quantum melting. How¬ 
ever, matter degrees of freedom at integer filling factors in the quantum Hall 
regime tend to be often well described by mean-field theory [37]. The accu¬ 
racy of mean-field theory is generally related to the restricted Hilbert space 
of LLs, which preclude the formation of competing correlated states with 
larger quantum fluctuations. There are several examples of interesting bro- 
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ken symmetry states in both semiconductor quantum wells and graphene that 
are accurately described by mean-field theory, including spin-polarized ferro¬ 
magnetic states at odd filling factors |60], coherent quantum Hall bilayers in 
semiconductors systems with coupled quantum wells |61j . and spin-density 
wave states in neutral graphene [3S]- In some cases, the state selected by 
mean-field-theory energy minimization is the only state in the quantum Hall 
Hilbert space with a given set of quantum numbers, and therefore is exact. 
The situation here is similar to the coherent bilayer state [61] in that we have 
coherence between adjacent LLs. 

Finally, we mention that physics similar to that described in this Article 
is not expected to be limited to graphene but should equally occur in 2D 
electron gases in semiconductor (e.g. GaAs) quantum wells. There are a 
number of quantitative differences in detail, however. Most critically, the 
anharmonic LL spectrum of graphene should make it possible to achieve a 
better selective coupling to a particular M, M + 1 doublet [25] . 

Data availability. The data files used to prepare the figures shown in 
the manuscript are available from the corresponding author upon request.. 
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(a) 





Figure 1: Phases of a graphene integer quantum Hall polariton fluid. 

Pictorial representation of the two phases supported by a graphene quantum 
Hall fluid interacting with a uniform electromagnetic held (in yellow). Panel 
(a) When electron-electron interactions are weak, the ground state of the sys¬ 
tem l"^) is a spatially-uniform polariton condensate. In pseudospin magnetic 
language, this state is a collinear ferromagnet with all pseudospins, deflned 
in Eq. ([^, denoted by blue arrows in the Bloch sphere, pointing along a com¬ 
mon direction. Panel (b) When electron-electron interactions are sufficiently 
strong, the ground state of the system \'ip) spontaneously break translational 
invariance. In pseudospin magnetic language, this state is a spiral pseudospin 
state. The spiral is driven by antisymmetric interactions between z and x 
pseudospin components as explained in the main text. 


22 


(a) 


(c) 


e 


1.8 

1.6 

- ] —1 - 1 - 3 - i - 1 - 

1 /N / 

1.6 

- pn - 1 - 1 -^^- 

■ / 

1.4 

- 1 > y 

1.4 

- 1 

1.2 

- 1 1 

1.2 

■ / / 

1.0 

-I ' v -- 

X 1-0 

- 1 ' 

0.8 

1 ■ 

n 

o 

bo 

- 1 ' — "---^ 

0.6 

- 1 1 ' / 

0.6 

- / 

0.4 

-1 1 ' .■ 

0.4 


0.2 

. / / ■ 

• V 

0.2 

-1 ' ^ 

0.0 

f * — " 1 — — — 1 — — - 

0.0 

1 - 1 _ 1 _ 1 _ 1 _ 



5 0 5 10 15 20 25 

A/ff 

5 0 5 10 15 20 2 

A/5 

^ 1.6 

f 1 1 ' 1 1 1 

‘ 1 /■ 

1.6 

,1 ‘P ' 1 1 

~ 1 / 

1.4 

■ 1 / 

1.4 

- 1 / 

1.2 

- ' / . : 

1.2 


X 1-0 

- / • .*“ _ _: 

X 1-0 

- I • 

77 

o 

bo 

- 1 -W ." 

77 

o 

bo 

-1 ^ — — — — — 

0.6 

- 1 ; 

0.6 


0.4 


0.4 

.1 ;( 

0.2 


0.2 

; ' ; " _ 

0.0 

J _ 1 _ • 1* 1 — -p — — 1-^ — - 

0.0 

J _ £ _ Lj _ 1 _ ^^ _ 


10 


15 


20 


25 


-5 


10 

A/ff 


15 


20 


25 


Figure 2: Phase diagram of a graphene integer quantum Hall polari- 

ton fluid. White (grey-shaded) regions represent the values of the detuning 
A—in units of g —and density rix of total excitations at which the homo¬ 
geneous phase described by Eqs. (|T7|)-(@ is stable (unstable). Panel (a) 
Kr = 5 and M = 1; panel (b) = 5 and M = 2; panel (c) = 15 and 
M = 1; panel (d) Kj- = 15 and M = 2. In each panel, blue lines denote the 
location of points in the plane (A/^f, nx) where the ratio between the number 
of excitons and the number of photons is constant: rie^/npY^ = 1/10 (dashed 
line), rtex/rtph = 1 (dash-dotted line), and rtex/rtph = 10 (dotted line). These 
curves have been calculated with reference to the homogenous phase. 
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Figure 3: Collective excitation spectrum of the homogeneous fluid 
phase. Panel (a) Dispersion relations of upper and lower dressed polariton 
modes (solid lines) in the long-wavelength q^B ^ limit. The dashed line 
(dash-dotted line) represents the cavity photon dispersion (bright electronic 
collective mode), when the electron-photon coupling ^fq is set to zero. In this 
panel ipq = 0. The inset shows a zoom of the lower polariton dispersion 
relation for q^-Q —t 0. Note that the dispersion behaves as y/g in this limit. 
Panel (b) Dispersion relation of the lower polariton mode, along three direc¬ 
tions: (pq = (f) (solid line), ipq = (p + n/A (dash-dotted line), and pq = (f) + n/2 
(dashed line). All the data in this hgure have been obtained by setting = 5, 
B = 0.5 Tesla, M = 1, nx = 0.1, and A = g. 
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Figure 4; Variational parameters of the modulated phase. This figure 
shows the optimal values of the variational parameters u,v,6 and Q* in 
Eq. (30) for a cavity dielectric constant = 5, highest occupied LL M = 1 
and different values of the detuning A. Panels (a) and (c): A = 2g. Panels 

(b) and (d): A = 3g. Panels (a)-(b) Dependence of the characteristic wave 
number Q* (in units of l/is) on the density rix of total excitations. Panels 

(c) -(d) Dependence of the quantities u (blue squares), v (green triangles), 
and 6 (cyan diamonds) on the density nx of total excitations. Grey-shaded 
areas have the same meaning as in Fig. 
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Figure 5: Energetics and photon densities in the modulated phase. 

Data denoted by red circles represent the quantity Ve in Eq. (33) plotted as a 
function of nx- Data denoted by blue squares represent the ratio between the 
photon density Uph and the density rix of total excitations, as a function of 
nx and for the modulated phase only. Data in this figure have been obtained 
by setting the dielectric constant at Kr = 5 and the highest occupied LL at 
M = 1. Above the horizontal blue line, n^^/nx > 10%. Different panels 
refer to different values of the ratio A/^f; panel (a) A/g = 2, (b) 2.5, (c) 
2.75, and (d) 3. 
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Supplementary Figure 1: Dimensionless strength of electron-electron 
interactions. Dependence of Oee (in units of Oee^c) on the Landau level 
index M. 
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Supplementajy Figure_2: Pseudospin-pseudospin interactions. Fourier 
transforms Ji{q) and V 2 {q) (in units of Settee)- Panels (a), (b) and (c) refer 
to the symmetric interactions, while panel (d) refers to the antisymmetric 
interaction. The red solid (blue dashed) line refers to the case with highest 
occupied LL M = 1 (M = 2) in conduction band A = +. 
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Supplementary Figure 3: Phase diagram in the absence of symmetric 
interactions. The red dashed boundary of the grey-shaded area—which 
represents the unstable region for 57? 7 ^ 0 as in Figure 2(c) of the main text— 
enlarges: for = 0 (absence of symmetric interactions) the new boundary 
is represented by the blue filled circles. 
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Supplementary Note 1 

At low energies, charge carriers in graphene are modeled by the usual single¬ 
channel massless Dirac fermion (MDF) Hamiltonian [1, 2] 

■Hd = vt,(J ■ P , (1) 

where ud ~ 10® m/s is the Dirac velocity. Here a = {(Jx-,(yy) is a 2D vector 
of Pauli matrices acting on sublattice degrees-of-freedom and p = —ihVr 
is the 2D momentum measured from one of the two corners (valleys) of the 
Brillouin zone. 

A quantizing magnetic field B = Bz perpendicular to the graphene sheet 
is coupled to the electronic degrees-of-freedom by replacing the canonical 
momentum p in Eq. (1) with the kinetic momentum 11 = p -|- cAq/c, where 
Aq is the vector potential that describes the static magnetic field B. The 
corresponding Hamiltonian is 


'Ho = v-£>a ■ n . 


( 2 ) 


We work in the Landau gauge Aq = —By-k. In this gauge the canonical 
momentum along the x direction, px, coincides with the magnetic translation 
operator [3] along the same direction and it commutes with the Hamiltonian 
Hq. Thus, the eigenvalues oipx are good quantum numbers. A complete set 
of eigenfunctions of the Hamiltonian Hq in Eq. (2) is provided by the two 
component pseudospinors [4] 


(r|A,n, k) 


^ikx 


W_^n4'n-l{y - i%k) \ 

Xw+^n4’n{y - i%k) ) ’ 


( 3 ) 


where A = -|- (—) denotes conduction (valence) band levels, n G N is the 
Landau level (LL) index, and k is the eigenvalue of the magnetic translation 
operator in the x direction. In Eq. (3) 

w±^n = a/I ± 5n,0 (4) 


guarantees that the pseudospinor corresponding to the n = 0 LL has nonzero 
weight only on one sublattice. Eurthermore, 0n(|/) with n = 0,1, 2,... are the 
normalized eigenfunctions of a ID harmonic oscillator with frequency equal 
to the MDE cyclotron frequency Uc = v^ud/^b- Here = a/^c/ {eB) ~ 
25 nm/A/H|TbsIaJ is the magnetic length. 
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The spectrum of the Hamiltonian (2) has the well-known form [4] 

£x,n = XhUcx/n . (5) 

Each LL has a macroscopic degeneracy Af = NfS/{2TTi%) = where 

A^f = 4 is the spin-valley degeneracy and S' = is the sample area. 

The fully microscopic matter Hamiltonian is written as 


n 


mat 


Hq + He 


( 6 ) 


Here, Ko is the second-quantized version of Eq. (2), 

A,n,fc,^ 


(7) 


where (cx,n,k,i) is a fermionic creation (annihilation) operator for an 

electron with band index A, LL quantum number n, and eigenvalue of the 
magnetic translation operator along the x direction equal to k. The collective 
index ^ refers to the valley (iL, K') index and spin-projection along the z 
direction. The second term in Eq. (6), "Hee, represents Coloumb interactions. 
This term can be written as 


n 


ee 


1 


[n_qnq 

q 




( 8 ) 


where nq is the Eourier transform of the electronic density operator 


n, 


q = / dre , 


(9) 


and Uq is the 2D Eourier transform of the Coulomb potential 

27re^ 




eg 


( 10 ) 


Here Kj. is the cavity dielectric constant. The real-space field operators 
and V’(i') have been introduced in the main text. 

Selecting only terms that respect the SU(4) spin-valley symmetry, the e-e 
interaction term in Eq. (8) can be written in the following manner [4] 

«- = E k )V(Aini),(A2n2),(A3n3),(A4n4) (q) 

q ^i,ni i,^',k,k' 

^ ^Xl,ni,£^,k^\2,n2,i',k',k'^>'A,ni,i,k-q^ i (H) 


5 



where 


1^(Ai,ni),(A2,n2),(A3,n3),(A4,n4) (q) '?^q‘^(Ai,ni),(A4,n4) ( Q)‘^(A2,n2),(A3,n3) (q) • (12) 

The form factors ^(A,n),(A',n')(Q) given by [4] 


^(A,n),(A',n')(q) 


1 

2 


W——^n' Dn—l^n'—1 



W W+ ,n' DJi n' 



Here, q = Qx + iqy and q* = qx - iqy, 


Dn,n' (^) 



Dn',n{-Z*), 


for n>n' 
for n <n' 


(13) 


(14) 


and ^\x) are generalized Laguerre polynomials [5]. 

We consider the integer quantum Hall regime in which a given number 
of LLs are fully occupied and the remaining ones are empty. Since the MDF 
Hamiltonian is particle-hole symmetric, we can consider, without loss of gen¬ 
erality, the situation in which graphene is n-doped and the Fermi energy lies 
in the conduction band (A = -f). We denote by n = M the highest occupied 
LL, and the lowest empty LL is n = M -|- 1 We are interested in the case 
in which cavity photons with energy hw are nearly resonant with the energy 
difference between the two conduction-band LLs n = M, M + 1. In this limit, 
the fermionic Hilbert space can be reduced to the resonant doublet. From 
here on, we denote by c\^n,k,i and ^xnki operators for states which do not 
belong to the resonant doublet M,M + 1. We will keep using c\^n,k,i and 
c\nki states which belong to the resonant doublet. We can rewrite 

the full microscopic Hamiltonian in Eq. (6) as the sum of three terms: 


'Hmat = "Hd + "Hm + 'Hdm • 


(16) 


The first term "Hd contains only fermionic field operators related to the reso¬ 
nant doublet M, M+1. The second term, "Hm, contains only field operators of 
the type c\^n,k,i-, ^\nk(- these degrees of freedom play the role of a “medium” 
for the resonant doublet. The third term, "Hdin, describes coupling between 
the resonant doublet and the medium degrees of freedom. 
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To obtain an effective matter Hamiltonian, we start from the fully mi¬ 
croscopic Hamiltonian "Hmat and we treat in an exact fashion all the terms 
that involve field operators {c\^n,k,i and acting only on the doublet 

n = M, M -|- 1 in conduction band A = -f. All the other terms (containing 
c\,n,k,i and <^xnk^ treated within the Hartree-Fock approximation [3]. 
The medium Hamiltonian is discarded. In the coupling term "Hdin, we only 
keep terms that separately conserve the number of particles in the M, M + 1 
doublet and in the medium. These are terms of the form (fc^cc. Terms of 
the form c^c^cc are discarded. We therefore replace 

^Ai,ni,fciW,n2 Mfy W,n3 M 

'5Ai,A4<5ni,n4^fcl,A:4<5^l,.j4 0(^+,M ~ ^Ai,ni)c_|_ ‘ 

After straightforward algebraic manipulations we reach the final result 
for the effective matter Hamiltonian, which is best expressed in a pseudospin 
representation in which pseudospin “up” (“down”) corresponds to the M + 1 
(M) LL. To this end, we define the following pseudospin operators: 


rm 




n,n'=M,M+l k 


-il\qy{k^q^/2) 


(17) 


and 



(18) 


where m = 0,..., 3 and [Tm\nn' labels the matrix elements of a four-vector 
of 2 X 2 matrices acting on the M,M + 1 doublet, specifically Tq represents 
the 2x2 identity matrix and represent the ordinary 2x2 Pauli 

matrices. 

The final effective matter Hamiltonian reads as following: 


= El,V^S,{0) + ^ Knm'(q)^m(-q)^„i'(q) • 

_( 18 ) 

Here, Efj = F/m -I- Ej^q and = Hat T with Ej^ = huc (a/A f -|- IT -|- 

a/M)/ 2. The terms Em and are due to exchange interactions between 
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the resonant doublet M, M + 1 and occupied LLs outside the doublet, i.e. 


Em — 


ZE' V(+,M+l),(A,n),(+,M+l),(A,n) (q) 


q A,n 


^(+,M),(A,n),(+,M),(A,n)l 


“ “972 V(+,M+l),(A,n),(+,M+l),(A,n)(q) 


q A,n 

- V(+,M),(A,n),(+,M),(A,n)(q) , (21) 

where the summation over A,n runs only values such that eA,n < 

The last term in Eq. (19) describes e-e interactions within the resonant 
doublet M, M + 1, i.e. 

HTtm'(q) = ^ 'y ) ['^■m]ni,n4 ['rm']n2,n3 V(+^ni),(+,n2),(+,n3),(+,n4) (q) • (22) 

ni=M,M+l 

This is a 4 X 4 Hermitian matrix, which can be decomposed into its real and 
imaginary parts: Knm'(q) = Re[h7r„'(q)] + *Im[V7m'(q)], where 


R.e[V7m'(q)] — R-G[hrn'm(q)] 
Im[l44m'(q)] = -lm[Kn'm(q)] 


It is possible to show that every quantity Knm'(q) is either purely real or 
purely imaginary: I4nm'(q) is purely imaginary if m = 1 or 2 (m = 0 or 3) 
and simultaneously m' = 0 or 3 (m! = 1 or 2); for any other value of m,m' 
7w(q) is purely real. 


Supplementary Note 2 

We consider a graphene sheet coupled to the electromagnetic field in a cavity. 
The Hamiltonian that describes the coupling between electrons and cavity 
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photons is 


+ ^ U^X'n'(^erai^y ^)^n',n-l\ {o,ci,u (i-q,u)'^X'n'k+q^^>^nk ) (24) 

where ^fq = hwc^/e^]~(^2eL^Ko^ is the light-matter interaction parameter, 
e^(q, z/) = (x ± iy) ■ eem(q, where eem(q, is a unit vector describing 
the linear polarization v of the electromagnetic field, cUq = 
is the cavity photon dispersion relation, is the cavity dielectric constant, 
V = LzL?' the volume of the cavity, and L is the cavity length along 

the z direction. 

Among all the processes described by "Hint, we take into account only reso¬ 
nant terms which describe the photon-induced electronic transitions between 
LLs M and M -|- 1 in conduction band. This approximation is called rotat¬ 
ing wave approximation (RWA). In the RWA, the light-matter interaction 
Hamiltonian becomes 


Ml„, = + “-q,LS'+(q) 


(26) 


where Oq^L = (flq,a; ” [ajj l = (®q,a; + *®q,y)/v^] the annihilation 

[creation] operator for a left-handed photon and ^±(q) = [*S'i(q) ±iS'2(q)]/2. 
The term containing 5*+ (S-) describes transitions from LL M (M + 1) to LL 
M + 1 (M) assisted by the annihilation (creation) of a left-handed photon. 
For B = Bz, the RWA selects the left circular polarization, eem(q, L) = 
(x — iy)/x/2\ of course, the RWA will select the right circular polarization, 
eem(q, R) = (x iy)/ v4, for B = -Bz. 

Within the RWA, the sum of the number of cavity photons iVph and the 
number of excitons Wx is a conserved quantity. The number of photons can 
be expressed in terms of the photon field operator as 

^Ph = > (26) 

q ,!2 

while the number the number of excitons can be expressed in terms of the 
pseudospin operators as 

7Vex=Ar[l + A3(0)/v(A(l/2 . (27) 
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Supplementary Note 3 


It is possible to show that the variational state \il)) introduced in Eq. (3) of 
the main text can be re-written as: 


where 


and 




p' = — \aa 
7 


0 L + j , 


k,^ 


Id = J-^5^1og(cos4/2) , 




(28) 

(29) 

(30) 

(31) 


7 = vW+W ■ (32) 

We note that the commutator between the operator dd and its hermitian 
conjugate d is given by: 


[d,d^] — (6*fc/2) j • 

(33) 

Calculating the expectation value of [d, d^ on the variational state introduced 
in Eq. (3) of the main text and taking the low-density limit 6*^ 1, we find: 


{ilj\[d, dd]\^jJ) ^ 1 


‘^Ekdl ■ 


(34) 


If the second term on the right-hand side of the previous equation is neglected, 
excitons can be treated as bosons. In the same limit, the operator pi has the 
meaning of a polariton creation operator. The variational state reported in 
Eq. (3) of the main text is a coherent state of polaritons and the quantity 
A/’lyp represents the average number of polaritons. 
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Supplementary Note 4 

We here report a number of relevant parameters that appear in the energy 
functional introduced in Eq. (7) of the main text: 

^ee = 'Yj} ^ W q I ^ ^ I ~ l‘^{+,Ar+l),(A,n) (q) I ] 

- -^(+,Ar),(+,Ar)(q)-^(+,M+i),(+,M+i)(q)| , (35) 



and 

q 

“ l‘^(+,M+l),(+,M+l) (q)P-|:f(+.M),(+.M)(q)l"] ■ (37) 

As in Supplementary Note 1, the sum over A,n runs only over indices such 
that eA,n < e+,Ar- We note that a^e involves only the two resonant LLs M 
and M + 1. 

Supplementary Figure 1 shows the quantity Oee (in units of Oee^c) as 
a function of the LL index M in the interval 1 < M < 15. Here a^e = 

e^/{K^hvY:)) ~ 2.2!Kr- 

The symmetric and antisymmetric pseudospin-pseudospin interactions 
can be decomposed as following: 


Jeik - k') = J^k - k') + Je,.ik - k') (38) 

and 

V2{k - k') = V2,dik - k') + V2,^{k - k') . (39) 
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The direct contributions are given by: 



>Jl,dik 

n 

(40a) 


J2,A{k 

H 

-k') = 0, 

(40b) 


J^A^k 

q 





(40c) 


'h>2,x{k 

q 




X [‘^(+,M+l),(+,M+l)(Q'y) -‘^(+,M),(+,M)(Q'y)] • 

(40d) 

The exchange ^ 

contributions are given by: 



-k') 

= 2^2 {-^(+,Ar+i),(+,Ar+i) (q) J'(+,m),(+,m) (q) 

n 




M 

- Re^ [j'(+,M+i),(+,M)(q)] , 

(41a) 

<J2,x{k 

-k') 

— 21,2 y^'^q{>^(+,Ar+l),(+,Ar+l)(q)-^(+,M),(+,M)(q) 

n 




M 

+ Re^ [j^(+^M+i),(+,M)(q)] }<5g^,fc-fc' , 

(41b) 


-k') 

^ 4 ^ y^'^q -^(+,M+l),(+,M)(q) — \J^i^j^^M+l),{+,M+l) 

q 

(q)r 




(41c) 

'h>2,x{k 

-k') 

— 21^2 (+,M+1),(+,M+1) (q) + (+,M),(+,M) (q)] 

n 




M 

X Re [j^(+^M+l),(+,Ar)(q)] , 

(41d) 

where Afs has 

been introduced in the main text. 


Supplementary Figures 2 (a), (b), (c) and (d) show the Fourier transforms 



^ /c 

(42) 

and 


-1 




^ /c 

(43) 
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for M = 1 (red solid line) and M = 2 (blue dashed line). In the g —)■ 0 limit 
we obtain: © 2 ( 0 ) = 0, JeiO) = j£,x(0) 7 ^ 0 for any i, and Ji(0) = ^ 72 ( 0 ). 
Moreover, we find that 

aee = 2[J3(0)-Jl(0)] . (44) 

Supplementary Note 5 

The correction Age in Eq. (35) to the cyclotron transition energy is related 
to the renormalization [6, 7, 8] of the Dirac velocity ud due to exchange 
interactions, which occurs also in the absence of a magnetic field. It is well- 
known [9, 10, 11] that Age is logarithmically divergent, i.e. 

Agg = ttgg^ [In(nmax) + Cm] , (45) 

O 

where Umax is a cut-off defined below and the constant Cm depends on the 
highest-occupied LL with index M. For example, Cq ~ —1.017 and Ci ~ 
—2.510. The Dirac model applies over a large but finite energy region, so we 
define a high-energy cut-off W in valence band. At any given magnetic field 
B, the integer Umax represents the number of LLs in the valence band with 
energy larger than W, i.e. 

nmax = ^ , (46) 

t> 

where Bw = eW"^/{2chvD), such that W = We follow 

Shizuya [11] in writing the correction to cyclotron transition energy in terms 
of the renormalized Dirac velocity: Dm = V^hv^/iB and = ud + ^ud- The 
quantity ud is the bare Dirac velocity and (Jud = c^qedc/{SK,^)[\ n{B]v/B) + 
Cm] is the correction to the Dirac velocity, where agso = e^/{hc) ~ 1/137 
is the QED fine-structure constant. This allows us to fix Bw = 450 Tesla 
to make sure that matches the value = 1T2 x 10® m/s measured [12] 
from the intraband 0 —)■ 1 cyclotron transition energy at = 18 Tesla and in 
a sample with = 5. We have fixed the bare Dirac velocity to v^, = c/300. 

Supplementary Note 6 

In this Supplementary Note we discuss the relative role of symmetric and 
antisymmetric pseudospin-pseudospin interactions in determining the precise 
form of the phase diagrams shown in Figure 2 of the main text. 
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We start by setting to zero all the symmetric interactions in the energy 
functional introduced in Eq. (7) of the main text: = 0. In this case, we 

find that the unstable regions (grey-shaded regions in Figure 2 of the main 
text) considerable expand. This is illustrated for the case = 15 and M = 1 
in Supplementary Figure 3. 


Supplementary Note 7 

In this Supplementary Note we explain more in detail the approach we have 
followed to find the elementary excitations of the polariton fluid. 

In the spin-chain language introduced in the main text, the state |V’) 
of the homogeneous fluid phase represents a collinear ferromagnet in which 
the expectation value of the spin operator S(q) = [*S'i(q), S' 2 (q), S' 3 (q)]"'" is 
non-zero only at q = 0 and oriented along the direction 

z' = [— cos(0) sin(6*), — sin(0) sin(6'), cos(6*)]'^ , (47) 

i.e. (V’IS(q)l'^) = —i'y/M The collection of units vectors x', y', z' with 

x' = [cos(0) cos(6'), sin(0) cos(6'), sin(6')]"'" , (48) 

and 

y = [- sin(0), cos(0), 0]^ , (49) 

forms an orthonormal set, which will be used below to construct low-energy 
collective excitations above the ground state |'^). 

To study the collective mode spectrum, we assume that \'4>) is subject 
to an infinitesimal pertubation, which induces a small change \5%Ij), i.e. we 
write 1'^') = 1'^) -|- The most general infinitesimal change orthog¬ 
onal to 1'^) can be written as a superposition of low-energy single-particle 
excitations [13] 

l<5^) = f5^«qWi,^ + 5^rq^P+(-q) J I'll)) , (50) 

V ?,q / 

where p^^(q) = p«^(q) ■ (x' + iy')/2 and p^^(q) = [pf (q), ^^(q), ^^(q)]^, 
with Pm(fl) in Eq. (17). In writing Eq. (50) we have assumed that the in¬ 
finitesimal pertubation does not induce any spin-flip or intervalley scattering 
process. 
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In the time-dependent Hartree-Fock (or generalized random phase) ap¬ 
proximation [13], the spin wave spectrum can be calculated by making use 
of the Heisenberg equation of motion (EOM) for the expectation value of an 
operator O, evaluated on the state |'^'), i.e. 

where the total Hamiltonian "H is reported in Eq. (1) of the main text. Writ¬ 
ing EOMs for O = P-(q), and p+(q) and keeping terms up to 
linear order in and = (pl^(q)), we find a homogeneous 

system of first-order linear differential equations. Eigenmodes are found by 
replacing ihdt —)■ fl. We obtain an eigenvalue problem, whose solution gives 
the spin wave spectrum. We find six independent collective modes, four 
due to the fourfold spin-valley degeneracy and the remaining two due to 
the two possible fight polarizations. Only two collective modes are hybrid 
in that they contain both light and matter components. They are com¬ 
posed by a left-handed photon OqL and a bright collective matter mode 
S'(_(q) = S(q) ■ (x' — iy')/2. In this subspace of hybrid modes, the eigenvalue 
problem reduces to (Oql — M)v = 0 where 1 is the 4x4 identity and the 
4x4 matrix M is reported in Eq. (22) of the main text. 
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